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We discuss, in quahtative and quantitative fashion, the yields of hadron 
resonances. We show that these yields, in general, are not in chemical 
equilibrium. We evaluate the non-equilibrium abundances in a dynamic 
model implementing the 1 + 2^3 resonance formation reactions. Due to 
the strength of these reactions, we show the I](1385) enhancement, and the 
A(1520) suppression explicitly. 

PACS numbers: 25.75.Nq, 12.38.Mh 

1. Why study strange resonances and quark gluon plasma? 

We study the quark confining vacuum structure, with the objective to 
liberate quarks and gluons at high temperature. The present day experi- 
mental effort involves colliding heavy nuclei and this limits the domain of 
deconfinement to the size of the atomic nucleus. The 'free' quarks and glu- 
ons form a thermal gas comprising color charges, the quark-gluon plasma 
— a hot soup of elementary matter, last seen in the Universe when matter 
was formed at about 30 fis after the big-bang. 

The temperature at which we deconfine the quark content of nucleons 
is Tc ~ 160 MeV. However, we reach in heavy ion collisions, Tyiax — 2-5 Tc. 
At this relatively high temperature, we create strange quark and antiquark 
pairs in an abundance which rivals that of light quarks [Ij. While this is 
going on, the compressed quark matter expands, the expansion consumes 
the thermal pressure and energy, and ultimately the fireball of hot quarks 
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breaks apart at Tj < Tc- Of particular interest is the high abundance of 
antiquarks, including anti-strangeness. The high strange antimatter yield is 
our evidence that we have recreated the early Universe in a laboratory set 
experiment. 

We are interested in understanding the physics of quark-gluon plasma in 
the last moments of its existence, and thus the conditions prevailing in the 
transformation of quarks and gluons into hadrons. The yields of particles 
produced can be successfully described using statistical physics methods. 
We use the program SHARE (Statistical HAdronization with REsonances) 
for this purpose [21 [3] . The remarkable success of this model indicates that 
the production of hadrons in heavy ion collisions is governed mostly by the 
accessible phase space. This, in turn, implies that particles are produced in 
a process resembling vapor evaporation from a hot soup. In this process, 
hadrons are formed from very 'sticky' quarks and this helps to saturate the 
probability of particle formation. 

Aside of stable (under strong interactions) particles, the statistical hadro- 
nization of quark-gluon plasma predicts the yields of (anti)baryon reso- 
nances. Experimental results available today show that the yields of these 
states do not always follow the model expectations [5l [H [3 [HI [9l [10] . Given 
the success of SHARE, we interpret this as a post hadronization evolution 
of observable yields [HI [12]. Thus, the first objective of this report, see 
section [2l is to explain why, for stable particles, we can directly apply the 
statistical hadronization model, while observed resonance yields are subject 
to post-hadronization dynamics. 

After that, we develop the kinetic model tools in section [3l define the 
model and the initial conditions in section [U and present, in the following 
section [5l a detailed numerical study demonstrating that S(1385) can be 
greatly enhanced in the observed abundance compared to statistical equilib- 
rium, while other more stable resonances, such as A(1520) are suppressed. 

The understanding of this behavior of strange baryons and antibaryons 
and their resonances sharpens the tools available to us in the study of quark- 
gluon plasma properties, at the time of phase conversion into hadrons. This 
work improves the understanding of the physics of hadronization of quark- 
gluon plasma, that is of the process of freezing of the deconfined vacuum. 

2. Why resonances, in general, do not chemically equilibrate? 

The observed stable particle yield is controlled solely by freeze-out tem- 
perature, and this yield contains the decay of all resonances. However, the 
resonance abundances can evolve and mix with stable particles without al- 
tering the observed final stable particle yield, since there is no information 
in the stable particle yield about 'sharing' of the yield with resonances. 
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Therefore, we cannot assume that resonance yields are governed by same 
physics as the final stable particle yields. We consider, as an example, the 
reaction 

A + vr^E*, (1) 

and imagine, for purpose of following simple illustration, that the system 
we study comprises ONLY these three particles. 

Reaction Eq. ([T|) does not change either the pion vr or lambda A yield, 
since all E*-resonances ever made will ultimately contribute to these yields 
upon resonance decay, which naturally happens before the stable particles 
are observed. On the other hand, we measure the yield of S* by the invariant 
mass method. This is done by considering all pair combination of two 
presumed decay products and evaluating from the energy and momentum 
of both particles the invariant mass distribution dN/dM where: 

m = ^{e^ + Ek?-{p^ + pk?. 

This method implies that as long as the decay products do not rescatter 
after decay, that is before leaving the medium, the yield of resonances is 
determined by the observed yield of the decay products. On the other 
hand, we can assume that each elastic scattering deflects the momentum 
vector, so that the invariant mass method fails to observe the resonance. 

Given this consideration, we can evaluate the relative yield of S*/A 
assuming that the system expands very slowly (and we ignore spin and 
isospin for simplicity). We are given the inelastic reaction rate, i?(T)in, 
at temperature T derived from the cross section governing the inelastic 
reaction Eq. ([1]). Similarly, we have total elastic rate, i?(T)ei, originating in 
any elastic scattering (that is of vr or A) in the medium. The elastic reaction 
rate is proportional to the probability of not observing a S*, while the 
inelastic reaction rate, with the same proportionality constant, is describing 
the probability of seeing a S*. This means that the invariant mass method 
observed relative S* yield is the ratio of inelastic scattering rate to any 
reaction rate in the medium, that is, 

_ R{Tr)\n 

A ~ i?(r,)in + i?(r,)ei' 

considering that the nucleon yield already comprises all S* ^ A decays. 
Here, we denote by the resonance free-streaming temperature which is 
typically lower than the stable particle freeze-out condition Tf. Note that 
if the elastic reactions are negligible, in this greatly simplified model all A 
are descendants of S* and the ratio is unity. 

We see, in this example, that the observed relative resonance yield is 
not governed by thermal/statistical properties of the medium, but depends 
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Fig. 1. Reactions scheme for A(1520) and S](1385) population evolution. 

decisively on the strength of the inelastic reaction rate derived from known 
resonance cross section of reactions, such as Eq. ([1]). In this work, we will be 
addressing realistic physical system and will focus our interest on the under- 
standing of the role of the inelastic reaction rate of important resonances. 
We will be following their abundance evolution as function of time t, that is 
for a given function, T{t), of temperature. We thus will present our results 
as a function of T and focus interest on the range Tr to Tj, where Tj is the 
hadron formation temperature in QGP breakup. 

3. Time evolution equations 

In figure [H we show the scheme of reactions which all have a noticeable 
effect on A(1520) yield after the chemical freeze-out kinetic phase. The 
format of this presentation is inspired by nuclear reactions schemes. On the 
vertical axis, the energy scale is shown in MeV. There are three classes of 
particle states, which we denote from left to right as W (S=0 baryon), 'S' 
(S = — 1, 1 = 1 hyperon) and 'A' {S = —1, 1 = hyperon). 

Near each particle bar, we state (on-line in blue) its mass, and/or angular 
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momentum, and/or total width in MeV. The states A(1520) and S(1385) 
are shown along with the location in energy of A(1520) + tt and S(1385) + 
TT respectively, both entries are connected by the curly bracket, and are 
highlighted (on-line in red). The inclusion of the vr-mass is helping to see 
the kinetic threshold energy of a reaction. The lines connecting the A'^, S 
and A columns are indicating the reactions we consider in the numerical 
computations. All reactions shown in figure [1] can go in both directions, 
as shown by the double arrows placed next to the numerical value of the 
partial decay width Fj, in MeV. 

The evolution, in time, of the resonance yield is described by the process 
of resonance formation in scattering, 1 + 2 3, less natural decay 3 1 + 2: 



1 dN^ 
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(3) 



where subscripts i, j denote different reactions channels when available. We 
further allow different sets of subscripts i, j in order to allow more complex 
dynamical cases in which not all production and/or decay channels are 
present. 

Allowing for Fermi-blocking and Bose enhancement in the final state, 
where by designation particles 1 and 3 are fermions (heavy baryons) and 
particle 2 is a boson (often light pion), we have for the two rates: 
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and in analogy, we have for the 3 back-production rate 
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where (5^(1 + 2 — 3) = 5'^{pi + p2 — Pz) assures 4- momentum conservation 
and Qi, i = 1, 2, 3 is particle degeneracy. The Bose function for particle 2, 
and Fermi distribution for particles 1, 3 are: 
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Here, Tj is particles fugacity, and u ■ pi = Ei, for u'^ = (1,0) in the rest 
frame of the heat bath where d'^pdo{pf — mf) d^pi/Ei for each particle. 
Hence, Eq.(j31) and Eq.Q are Lorentz invariant, and thus as presented these 
rates can be evaluated in any convenient frame of reference. Normally, this 
is the frame co- moving with the thermal volume element. 

Since particles 2, 3 are here heavy baryon (resonances), we can work 
using the expansion of the relativistic distribution, the first term is the 
Boltzmann limit: 

Y " '^i^dixfKiixi), (7) 

where Xi = rrii/T, K2{x) is Bessel function (not to be mixed up with parti- 
cle 2). However, we use the complete Bose distribution to describe pions. 
We introduce in medium lifespan of particle 3: 



1 _ Ei-Rl23 



(8) 



and, similarly, channel lifespan rj, omitting the sum The rate -R123 is: 



R is independent of the fugacity, in the Boltzmann- limit. 

The production and decay rates are connected to each other by the 
detailed balance relation [THl FH]: 

^1 ^2 = T3 = R123. (10) 

Using detailed balance Eq. (fTOll . we obtain for fugacity T3 the evolution 
equation: 



dr 
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where we have also introduced characteristic time constants of temperature 
T and entropy S evolution 

1 _ d\n{x^''K2{x;)) ^ 1 _ d\n{VT^) ^ ^^^^ 

tt dT ' Ts dT 

The entropy term is negligible, ts ^ T3, tt since we implement near con- 
servation of entropy during the expansion phase. We implement this in way 
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which would be exact for massless particles taking VT^ = Const.. Thus, 
there is some entropy growth in HG evolution we consider, but it is not 
significant. In order to evaluate the magnitude of tt, we use the relation 
between Bessel functions of order 1 and 2 (not to be mixed up with particles 
1, 2) d{z^K2iz)) /dz = z^Ki{z). We obtain 

1 Ki(x3) f 

— = -TT-r^x-i-, (13) 

Tt > 0. For a static system with tt 0, we see that Eq. (jlip has transient 
stable population points whenever 

^TlTyr|-T3^1/r^=0. (14) 



Finally, we consider the evolution in time of Ti and T2. In the equation 
for Ti, we have terms which compensate what is lost/gained in T3 see 
Eq. (jlip . Further, we have to allow that particle '1' itself plays the role 
of particle '3' (for example, this is clearly the case for A(1520)). That is 
accomplished introducing a chain of populations relations as follows: 

(l' + 2'^l) + 2^3. (15) 

In the present setting, T2=7r = Const.. By virtue of entropy conservation 
the same applies to the case 2' = tt. However, if either particle 2 or 2' is 
a kaon, we need to follow the equation for T2,2'=x which is analogous to 
equation for particle 1 or 1'. 

4. Model details and initial conditions 

The evolution equations can be integrated once the time dynamics of 
the fireball and the initial conditions are fixed: 

1) We choose a model of expansion which fixes the behavior T(t); here, we 
invoke a model of matter expansion where the longitudinal and transverse 
expansion is considered to be (nearly) independent: 



T_ lf2{vT/R±) + l 
T~~3 



(16) 



where R± is the transverse radius, v is the velocity of expansion in trans- 
verse dimension. All flow parameters (or temperature dependence on r) are 
the same as in [15l [16] . 

2) We determine the initial values of particle densities (fugacities) estab- 
lished at hadronization/chemical freeze-out. We determine these for RHIC 
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head-on Au-Au collisions at y^SNN = 200 GeV. We introduce the initial 
hadron yields inspired by a picture of a rapid hadronization of QGP in which 
quarks combine into final state hadrons. For simplicity, we assume that the 
net baryon yield at central rapidity is negligible. Thus, the baryon-chemical 
and strangeness potentials vanish. The initial yields of mesons {qq, sq and 
baryons qqq, qqs are controlled aside of the ambient temperature T, by the 
constituent light quark fugacity 7g and the strange quark fugacity 7^. 

3) Since high energy pions are moving faster than the bulk of the matter 
and leave the domain in which the slower baryons are found, we assume that 
it is impossible to excite reactions with high threshold energy. We thus ex- 
clude channels for resonance 3 production with threshold energy AE > 300 
MeV. 

4) We characterize the hadronization dynamics: we assume that the strange- 
ness pair-yield in QGP is maintained in transition to HG. This fixes the 
initial value of 7^. In fact, since we investigate relative chemical equilibrium 
reactions, our results do not depend significantly on the exact initial value 
■js and/or strangeness content. The entropy conservation at hadronization 
fixes 7q. For hadronization temperature T(t = 0) = Tq = 180 MeV, 7^ = 1. 
However, when Tq < 180 MeV, 7^ > 1 in order to have entropy conserved at 
chemical freeze-out. At Tq = 140 MeV, jq = 1.6 that is close to maximum 
possible value of 7^, defined by Bose-Einstein condensation condition [T7]. 

5) The initial particle yields are fixed in terms of fugacities: 

T?i=y)=7^s, T02=.)=7^ (17) 

or 

^1l=N)=ll'^l2=K)=llls, (18) 

where y = E, A is a hyperon, the particle 1 is a baryon and particle 2 is a 
meson. The particle 3 is always a strange baryon: 

Tl,^y) = 7^s. (19) 

Note that for 7^ > 1, we have always initially 



TiT 



= 7' > 1 • (20) 

t=o 



As a consequence, initially the pair of particles 1, 2 reacts into 3. 
6) We do not need to follow the evolution in time for the pion yield, which 
is fixed by conservation of entropy per unit rapidity, as incorporated in 
Eq. (|16p . Thus, it is (approximately) a constant of motion. This can be 
seen recalling that the entropy per pion is nearly 4 within the domain of 
temperatures considered. The conservation of entropy implies that pion 
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number is conserved, and with VT^ ~ Const.. This further impUes that, 
during the expansion, 

= 7g = Const., 
which we keep at the initial value. 



5. S(1385) and A(1520) yield results 

In figure [21 we present the fractional yields i;(1385)/Atot (left), and 
A(1520) /Atot (right) as a function of temperature of final kinetic freeze-out 
T. The results for the hadronization temperatures Tq = 140 (blue lines). 
To = 160 (green lines) and Tq = 180 MeV (red lines) are shown. 



0.5| r 




Fig. 2. The ratio I](1385)/Atot, on left, and A(1520)/Atot, on right, as a functions 
of temperature T{t) for different initial hadronization temperatures Tq = 140, 160 
and 180 MeV (blue, green and red lines, respectively, recognized also by their initial 
value along the hadronization curve (dot-dashed). 



In figure [51 the green dash-dotted line is the result when the kinetic 
freeze-out temperature T coincides with the hadronization temperature Tq. 
There is no kinetic resonance evolution phase in this case, only resonances 
decay after hadronization. This result is similar to SHARE result (purple, 
dotted line). The small difference is mainly due to us taking into account 
the decays 

S(1670, 1750) ^ A(1520) + vr, (21) 
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which are expected/predicted in [18j. Similarly, for S(1385) our results for 
Tq = T are different from SHARE results because we include the decay 

S(1670) ^ S(1385) + 7r, (22) 

which is part of recently updated particle data set [TU] . 

For all initial hadronization temperatures, as the freeze-out tempera- 
ture decreases, the suppression for A(1520)ob/Atot ratio is larger than for 
A(1520)/A(1520)o (at the same temperature T of final kinetic freeze-out). 
The effect is due in part to the asymmetry to excite S(1775) and heavier 
E* directly due to its high mass. The resonance yield suppression effect is 
approximately of the same magnitude for all hadronization temperatures 
Tq. However, the initial hadronization yield of A(1520) is sensitive to tem- 
perature, and decreases rapidly with T. Therefore, only for Tq = 140 MeV, 
a kinetic freeze-out temperatures « 95-105 MeV the ratio Aob(1520)/Atot 
reaches the experimental domain Aob( 1520)/ Atot < 0.042ib0.01 [5l[6] shown, 
in figure [2l by dashed lines. For the same initial conditions, that is for 
To = 140 MeV, we find [15, il6j the ratio S(1385)/Atot ~ 0.45 at T ^ 100 
MeV (and for the entire range 95-135 MeV, in good agreement with exper- 
imental data [11 [7] ) . 

6. Conclusions 

We find that the resonant hadron states, considering their very large de- 
cay and reaction rates, can often interact beyond the chemical and thermal 
freeze-out of stable particles. Thus, the observed yield of resonances is fixed 
by the physical conditions prevailing at a later breakup of the fireball mat- 
ter rather than the production of non-resonantly interacting hadrons. This 
study quantifies the expectation that, in a dense hadron medium, narrow 
resonances are 'quenched' |llj. 

Despite a scenario dependent resonance formation or suppression, the 
stable particle yields used in the study of chemical freeze-out remain almost 
unchanged, since all resonances ultimately decay into the lowest 'stable' 
hadron. Therefore, after a description, e.g., within a statistical hadroniza- 
tion model of the yields of stable hadrons, the understanding of resonance 
yields is a second, and separate task which helps to establish the consistency 
of our physical understanding of the hadron production process. 

Our results show that the observable ratio A(1520)ob/Atot can be sup- 
pressed by two effects. First A(1520) yield is suppressed due to excitation of 
heavy 5]*s in the resonance scattering process. Moreover, the final A(1520)ob 
yield is suppressed, because S*s, which decay to A(1520), are suppressed 
at the end of the kinetic phase evolution by their (asymmetric) decays to 
lower mass hadrons. 
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We resolve resonance puzzle in that we find that some resonances can be 
enhanced and some suppressed. Specifically $](1385) is strongly enhanced, 
since the dense pion gas especially for 7^ > 1 pushes the A into S(1385). On 
the other hand, the narrow A(1520) is depleted by pions pushing it over to 
high mass resonances, which later can decay without repopulating A(1520). 
This effect is particularly strong if we observe that there are fewer high en- 
ergy particles than Boltzmann distribution predicts in a rapidly expanding 
and cooling fireball. 
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